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Theory of scanning tunneling spectroscopy of low energy quasiparticle (QP) states in vortex lattices of d- 
wave superconductors is developed taking account of the effects caused by an extremely large extension of QP 
wavefunctions in the nodal directions and the band structure in the QP spectrum. The oscillatory structures in 
STM spectra, which correspond to van Hove singularities are analysed. Theoretical calculations carried out for 
finite temperatures and scattering rates are compared with recent experimental data for high-T c cuprates. 



The electronic structure of the mixed state 
in d-wave superconductors reveals a number of 
fundamentally new features (see jl]-[| and ref- 
erences therein) as compared to the case of s- 
wave compounds, where low lying quasiparticle 
(QP) states are bound to the vortex core and 
are weekly perturbed by the presence of neibour- 
ing vortices at magnetic fields H -C H C 2- The 
vanishing pair potential in the nodal directions 
results in the extremely large extension of QP 
wavefunctions which are sensitive to the super- 
fluid velocity (V s ) fields of all vortices and, thus, 
the electronic structure is influenced by the vor- 
tex lattice geometry. The resulting peculiarities 
of the local density of states (DOS) can be de- 
tected, e.g., by a scanning tunneling microscope 
(STM). In this paper we focus on the theory of 
scanning tunneling spectroscopy of low energy 
QP states in vortex lattices of d-wave supercon- 
ductors and compare the theoretical calculations 
with recent experimental data |?],|^] for high-T c 
cuprates, where the dominating order parameter 
is believed to be of d-wave symmetry. Hereafter 
we assume the Fermi surface (FS) to be two- 
dimensional (2D), take the gap function in the 
form Ak = 2Aok x k y /kp (the x axis makes an an- 
gle 7r/4 with the a axis of the C11O2 planes). Let 
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us orient H along the c axis (H cl <C H <C H c2 ) 
and consider two types of vortex lattices: (I) rect- 
angular lattice with primitive translations ai = 
axo, Sl2 = cayo! (II) centered rectangular lat- 
tice with ai = axo, &2 = »(xo/2 — cxyo), where 
Haa 2 = 4>o is the flux quantum, and xo, yo, zo 
are the unit vectors of the coordinate system. 

Van Hove singularities. Our consideration 
is based on the analysis of the Bogolubov-de 
Gennes (BdG) equations for low energy excita- 
tions with momenta close to a certain gap node 
direction (e.g., k x = k F x ): (H + H')g = eg, 
where g = (it, v) is the QP wavefunction, 
Hq = V F a z p x + VA&xPy, o x ,b z are the Pauli ma- 
trices, H' = MV F V SX (1 + a z ) + MV A V sy a x , M is 
the electron effective mass, p = —ihV — eA/c, 
V F = hkp/M, V A = 2A /{hk F ), H = -iJz , 
A = Hyx , V s = (V sx ,V sy ). The spectrum of 
the Dirac Hamiltonian H can be obtained using 
the usual quantization rule for a cyclotron orbit 
(CO) area HQ. The periodic potential H' re- 
moves the degeneracy of the discrete energy levels 
with respect to the CO center and induces a band 
structure in the spectrum . The general so- 

lution can be written in the form of a magnetic 
Bloch wave: 

g = Y^ e^qx+ZimM+Zinaqyafyy _ 2n oa, q), (1) 
n 

where n is an integer, and q is the quasi- 



momentum lying within the first magnetic 
Brillouin zone (MBZ): -ir/(2a) < q x < n/(2a), 
— 7r/(2<ja) < q y < Tr/(2aa). The wavefunction 
G(y, q) is localized in the domain with the size 
L determined by q and energy values. The po- 
tential H' results in the splitting of the CO near 
MBZ boundaries (see Fig. [l]) and the spectrum 
consists of branches which correspond to the split- 
ted portions of the CO. For large Dirac cone 
anisotropy q: = Vf/Va^ > 1 (a = kp^o/2) and 
e < 0.5e* (e* = nhV F /a ~ A 0y /H/H c2 ) the har- 
monics in Eq. (jl|) do not overlap (L < 2aa) 
and one can replace H' by the potential (H') x 
averaged in the x direction (see Such 
a simplification is a natural consequence of a 
small size of the cyclotron orbit (COl in Fig. |l|) 
in the nodal direction as compared to the size 
of the MBZ. The energy spectrum consists of 
branches e n {q x = irQ/a) = e*E n (Q 1 aa), which 
are displayed in Fig. g in the first MBZ for 
iraa = 50 and iraa = 100. The number of en- 
ergy branches which cross the Fermi level can be 
determined as follows: N ~ 2qy8q\\ ~ 2y/iraa, 
where gjj is the minimum possible size of the 
CO in the q\\ direction and 8q\\ is the distance 
between MBZ boundaries. Each energy branch 
has an extremum as a function of the momen- 
tum q x near the MBZ boundary at a certain e n 
(we neglect here additional extrema which ap- 
pear due to the exponentially small splitting of 
energy levels near the points of intersection of 
the branches in the E — Q plane). Due to the 



one-dimensional (ID) nature of the low energy 
spectrum the divergent contributions to the DOS 
take the form: SN(e) ~ \e — e n \^ x ^ 2 (e > e„ for 
energy minima and e < e n for maxima). The 
distance between these peaks 5e ~ e* /(2aa) co- 
incides with a characteristic energy scale corre- 
sponding to van Hove singularities which occur 
when the CO intersects MBZ boundaries in the 
q y direction (see Fig. ||). The crossover between 
ID and 2D regimes in the band spectrum occurs 
at e c ~ 0.5e*, when the CO size in the q± direc- 
tion becomes larger than the size of the first MBZ 
(C02 in Fig. [l]). For e > e c the q^-dependence of 
energy becomes essential and results in the ap- 
pearance of 2D critical points, i.e. 2D local max- 
ima (or minima) and saddle-points. Thus, instead 
of square-root van Hove singularities we obtain 
a set of discontinuities and logarithmic peculiar- 
ities (8N(e) ln\l - e/e n \) in the DOS, re- 
spectively. Obviously, these 2D singularities are 
more sensitive to temperature and finite lifetime 
effects and, consequently, the suppression of the 
corresponding oscillatory structure in the DOS 
should be stronger in the high energy regime. The 
above analysis can be generalized for gap nodes 
at k = ±kpyo: the corresponding energy scales 
take the form 5e ~ 0.5e*/a, e c ~ 0.5e*/<r. 

Even in the low energy regime the DOS oscilla- 
tions with the energy scale 5e are surely smeared 
due to a finite scattering rate T and temperature 
and can be observed only for a moderate Dirac 
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Figure 1. Cyclotron orbits (COl, C02) and MBZ 
boundaries for a square lattice, (qu, q±) defines a 

coordinate system whose origin is at the node, Figure 2. Energy branches for iraa = 50 (a) and 

with q± (q\i) normal (tangential) to the FS. naa — 100 (b) 



cone anisotropy and rather large magnetic fields. 
Comparing our results with a numerical solution 
H of the BdG equations for a = 1 and a = 5/2 
we find that the above mechanism gives a good 
estimate of the energy scale of the double-peak 
structure in the tunneling conductance at the core 
center at H/H c2 — 0.3 (Se ~ 0.1 A ) and can ex- 
plain the absence of this structure at low fields 
H/H c2 = 0.05 due to temperature broadening 
(T = 0.1T C > Se ~ 0.05A ). In principle, van 
Hove singularities may account for peaks with a 
large energy gap ~ Ao /4 observed experimentally 
at the vortex centers in YBaCuO at H ~ 6T 
provided we assume a ~ 1. Unfortunately the 
latter assumption is not consistent with the re- 
sults of thermal conductivity measurements || 
(a ~ 14), and, thus, the nature of the experi- 
mentally observed peaks is still unclear. It is also 
necessary to stress here that the critical points 
in the DOS are a direct consequence of perfect 
periodicity and the introduction of rather strong 
disorder surely remove these singularities. 

Zero-bias conductance. Hereafter we neglect 
the DOS oscillations, discussed above, and con- 
sider the peculiarities of the zero-bias tunneling 
conductance g(r) starting from a modified semi- 
classical model proposed in Q . According to this 
approach, the Doppler shift of the QP energy, 
which plays important role for e S Ao^/ H/ H C 2, 
appears to be averaged in the nodal direction due 
to an extremely large size of a semiclassical wave 
packet in this energy interval. Within such an ap- 
proximation a diagonal (retarded) Green's func- 
tion can be written in the form: 



G*(k,e,r) 



e + iT - hk F V a 



(e + iT-hk F V av ) 2 



,(2) 



where et is the normal state electron dis- 
persion, V av = (V s ) x + (V s ) y . The 
scattering rate T should be determined self- 
consistently: T — N(T, e)/(2N F r) (Born limit), 
T = N F T U /N(T, e) (unitary limit), where 2r and 
N F are the relaxation time and DOS at the Fermi 
level in the normal state, T u = ni mp / (irN F ), Ui mp 
is the concentration of of point potential scatter- 
ed, and N(T,e) = -Im J G R d 2 k/(2ir 3 ) is the lo- 
cal DOS. Let us first consider the effect of finite 
temperature on the zero-bias conductance in the 



clean limit QT — ► 0). The expression for the nor- 
malized conductance reads: 



<?( r ) 



Tlnl 
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gN 



-Hog 



N(T = 0, e)de 
4:N F Tcosh 2 
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4T 



(3) 



A 2A ^ 
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Here gN is the normal state conduc- 
tance, $ x = <£>(x/R x ), $ y = $(y/Ry), 
$(z) = 2z — (2m + 1) for m < z < m + 1 (m is 
an integer), R y (R x ) is the distance between 
the lines parallel to the x (y) axis and passing 
through the vortex centers, e* = irhV F HR x /(f>o, 
e* = nhV F a/R y , For type I (II) lattices we have 
R x — a, R y = aa (R x — a/2, R y — aa). One 
can separate two qualitatively different regimes 
in the behavior of the conductance: 
(i) superflow dominated regime T <C e* x y , 



(ii) temperature dominated regime T » e* 9 
. Tln2 irA aH 



fJ 



32TH c2 



(4) 
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$ y | and 



where F x {x,y) = \<S> X \{R X /R y ) 
F 2 (x,y) = <P 2 X (R X /R y ) 2 + i n F ig. | we dis- 
play the contour plots of the functions F\{x,y) 
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Figure 3. Contour plots of the functions F\ 
(a) and F% (b) which determine the spatial varia- 
tion of the zero-bias tunneling conductance for a 
square lattice (x' — x/a, y' = y/a). 



F2(x,y) for a square lattice of type I (which 
is close to the one observed experimentally in 
YBaCuO 0). There are two consequences of an 
increase in temperature: (i) first, the spatial di- 
mensions of peaks in the local DOS become rather 
small comparing to the intervortex distance only 
for T > T* ~ A Qy / H/H c2 ; (ii) second, the ampli- 
tude of the peaks appears to be essentially sup- 
pressed in the limit T ^> T* . For magnetic fields 
H ~ 6T (which is typically the field of STM ex- 
periment @J|) one obtains T* ~ 2QK. Thus, we 
conclude that the finite temperature effects can 
not explain neither the narrow zero-bias conduc- 
tance peaks observed in YBaCuO nor the absence 
of these peaks in BiSrCaCuO at T = A.2K. To 
explain these experimental facts it is necessary 
to take account of the finite lifetime effects which 
can stronly influence on the behavior of the DOS, 
as it follows from the results of Refs. |lO|-|li| ob- 
tained on the basis of the usual semiclassical ap- 
proach with a local Doppler shift. Starting from 
the modified semiclassical model (|J) we obtain 
the following expression for the tunneling conduc- 
tance at T = 0: 



iv(r,£ = o) 
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Obviously Born scatterers result only in a mod- 
erate change of the DOS (see [[To)) since the 
corresponding T value for Aqt 3> 1 is very 
small comparing to e* and the conductance is 
given by Eq. On the contrary, in the uni- 

tary limit the expression (|J) is valid only in the 
clean case T u -C T* y ~ O.le^ /Ao (for a square 
lattice T* x y ~ 0.1A H/H c2 ). In the dirty limit 
r« 3> r* we obtain: 



g~g(H = 0)ll 



64T UJ ff c2 



F2(x 1 y) 



(7) 



where g(H = 0) ~ 0.5^/T u / Aq. In the vicinity 
of each vortex center the local DOS exhibits a 
fourfold symmetry with maxima along the nodal 
directions in a good agreement with numerical 



calculations based on the Eilenberger theory [fL3| . 
For H — 6T finite lifetime effects become sub- 
stantial if we assume T u > 10~ 2 Ao. Thus, our 
approach allows to explain rather narrow con- 
ductance peaks (see Fig. 0b) observed near vor- 
tex centers in YBaCuO Q], even without taking 
account of the nontrivial structure of the tun- 
neling matrix element, discussed in p2[ . With 
a further increase of the T u value the ampli- 
tude of the peaks at the vortex centers vanishes: 
Sg ~ yj A /T U (H/H c2 ). Such a high sensitivity 
of the 5g value to finite lifetime effects can prob- 
ably explain the difficulties in the observation of 
these peaks in the mixed state of BiSrCaCuO ||| . 
Note in conclusion that according to Eq. (0) the 
spatially averaged DOS in the dirty limit varies as 
H rather than HlnH (the latter dependence has 
been predicted in within the semiclassical 

approach taking account of the local V s value). 
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